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We investigate the formation of an infinite cluster of entangled threads in a (2+1)- 
dimensional system. We demonstrate that topological percolation belongs to the univer- 
sality class of the standard 2D bond percolation. We compute the topological percolation 
threshold and the critical exponents of topological phase transition. Our numerical check 
confirms well obtained analytical results. 



I. INTRODUCTION 



In this paper we would like to link together two different subjects, percolation and topology, 
which since now were considered separately in physical literature. The problem discussed below 
deals with entanglement of threads: we apply to the threads a random set of operators of "wind- 
ings" and check whether two initial sets of threads become mutually entangled. The transition 
between entangled and nonentangled states may be considered as a "topological phase transition". 
We investigate the critical properties of such transition and show that it can be mapped onto 
the percolation problem. We also numerically investigate the critical exponent of the topological 
percolation and demonstrate that it is equal to the one of the two dimensional percolation. 

Let us remind the general formulation of the bond percolation problem. Consider a planar 
(for example, square) lattice of size L x L (L — >■ oo) with elementary cell of size 1x1. Each 
bond on the lattice is occupied with the probability p independently on all other bonds (i.e. with 
the probability 1 — p the bond is left empty). There exists a critical value p = pc sX which the 
infinitely large connected cluster of occupied bonds is formed. In other words, above pc at least 
two bonds belonging to opposite sides of the lattice are connected by a path of occupied bonds. 
It is known that for the two dimensional bond percolation Pc = \ The 2D percolation has 
been extensively studied during last decade. The scaling behavior of crossing probability (the 
probability to percolate in, say, horizontal direction as a function of p) in the vicinity of pc is 
defined by the correlation length exponent jH, A deep insight of scaling at the critical point 

has been achieved using the conformal approaches 

0,110,110,13 

The problem discussed below deals with entanglement of threads of finite length growing from 
a two-dimensional substrate. Apparently for the first time the continuous version of the same 
problem has been treated in Let us specify the model under consideration. Consider the two- 
dimensional square lattice of size L x L, serving as the "substrate" from which the threads start to 
grow. Each time moment one pair of nearest neighboring (along x- or y-axes) lines on the lattice 
may produce the full-turn entanglement. We assign the "generators" gx, g^^ to the full clockwise 
and counterclockwise turns along x-axis, and gy, g~^ to the full clockwise and counterclockwise 
turns along y-axis, and as it is schematically shown in the figure ^ 

Let us choose randomly any pair of nearest neighboring threads on the lattice and entangle 
them with the probability q; with the probability 1 — q we leave the pair of threads unentangled. 
If we decide to entangle the threads, then we do that clockwise with the probability r"*" and coun- 
terclockwise with the probability = 1 — Later on we shall consider only two extremal cases: 
i) r"*" = 1, meaning that all neighboring threads are entangled clockwise overall the lattice, and ii) 
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Figure 1: Generators of elementary entanglements of neighboring threads. 



r"*" = 2, denoting the situation when each entanglement may be either clockwise or counterclock- 
wise with equal probabilities. Let us stress that for any 7^ 1 there is the possibility to unwind 
the neighboring threads, what essentially complicates the problem. The general structure of the 
bunch of entangled directed threads is shown in the FigE) 
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Figure 2: Bunch of entangled threads. 



Now we are in position to formulate the main question of our interest. Suppose that each time 
step we apply randomly one clockwise or counterclockwise winding to the top ends of randomly 
chosen pair of neighboring threads on the lattice of size L x L. When time passes, more and 
more threads become topologically connected (entangled). We would like to compute the typical 
(critical) number of applied windings at which the threads on the opposite sides of the substrate 
form the topologically connected cluster, i.e. become entangled. This model sets the concept of 
"topological percolation". 



II. ENTANGLEMENTS AND THE LOCALLY FREE GROUP 



The generators Qx, 9x1 9y^ 9y turns are defined for each pair of nearest neighboring 

threads and, hence, should depend on the lattice coordinate (in FigQ the corresponding lattice 
indices are omitted for simplicity): Qxihj) ^ a clockwise (+)-generator and g~^{i,j) - a counter- 
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clockwise (— )~generator of threads and {i + along x-axis; gy{i,j) - a clockwise (+)- 

generator and gy^{i,j) - a counterclockwise (— )-generator of threads and (i + along 

y-axis. The pair of indices (i, j) is attributed to the {x, y) coordinates of a left or a bottom threads 
in a pair; i and j run along x and y axes correspondingly. 

The whole set of generators {gx{i,j), ^{hj), 9yihj), dy^i'Ll)} (for all 1 < {i,j} < N - 1) 
sets the surface locally free group LT^d [l^i ll^ • It is known [lj| that each generator g can be 
represented in a form g = (we have omitted indices for brevity), where o is the braid group 
generator. The generators g obey the following commutation relations: 

f g%{i,j)g%(i: = g%{^ ,j')g%{i,j) if > 1 or \j - j'\ > 

9yiiJ)9yii',f) = 9yii',j')9yihj) if K " I > or \j -j'\ > 1 

J)5^'(^',/) = 9^'{i'J')9Ai,j) if i - i' / 0, 1 or / - i / 0, 1 

. 9'x{i,j)9x^{'i',j) = 9l{hj)9y'{hj) = 1 

where we have attributed s = " + " for g and s = " — " for g^^ . 

We can reformulate our geometrical problem of entanglement of threads in terms of random 
walk on the group CJ-2D- Consider the square lattice of size L x L. There are = 2L^ threads 
passing through the vertices of the lattice along z-direction normal to the plane (xy). Here and 
later we assume the square geometry, but we can consider any other lattices in the same way. 
We generate the random word in terms of "letters" -generators of the group CJ^2D- Namely, we 
randomly apply one after another generators of CJ-2D (with appropriate probabilities) to the open 
ends of threads, respecting the commutation relations Let us repeat that the generators acting 
along X- and y-axes have equal probabilities. The (+) and (— )-generators have correspondingly 
the probabilities r+ and r~ = 1 — r+. 

When a random sequence of M "letters" is generated, we check whether there are at least two 
threads on opposite sides of the lattice which are topologically connected to each other (i.e. whether 
they belong to the same cluster of entangled lines). We are interested in statistics of such cross- 
lattice entanglements averaged over different realizations of random sequences of "letters". As one 
sees later, this problem has straightforward relation to the two-dimensional bond percolation on 
the same lattice. 



III. GEOMETRY OF ENTANGLEMENTS AND BOND PERCOLATION 



Let us describe a simple geometric model which visualizes the commutation relations of the 
surface locally free group, making them very transparent. We associate the "white" cells to the 
generators gx{hj) and gy{i,j) and "black" cells to the inverse generators g~^{i,j) and gy^{i,j)- 
We drop cells randomly, one cell per each time step. The commutation relations specify the 
geometry of the heap constituted by falling cells. For example, if the falling cell and one of the 
previous cells have only one common edge, the upper cell remains on the lower one; if however the 
cell falls strictly on the cell with opposite color, these two cells annihilate. The figure Q shows 
the particular example of two pairs of non-commuting (left) and commuting (right) generators. 

The general view of growing heap of white-black cells is shown in Figll] This picture demon- 
strates the typical configuration of cells for the system size L = 10 and = 0.5. 
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Figure 3: Visualization of some commutation relations in the group CJ-2D- 




Figure 4: Typical heap of colored cells representing clusters of entangled threads. 

Remind that any white (black) cell corresponds to the full clockwise (counterclockwise) turn of 
neighboring vertical threads. The empty horizontal intervals between cells mean that the threads 
in these intervals are not entangled. 

Consider now the projection of the heap of cells onto the horizontal xy-plane (the base surface). 
In this projection the single thread is represented by a dot (a lattice site). Each falling cell corre- 
sponds to the bond between two neighboring sites. If at least one cell in a projection (independent 
on the cell color) covers some bond, we close this bond. Otherwise we leave this bond open. 

It is not difficult to establish the bijection between entanglements and the standard percolation. 
Two distant threads are linked together (entangled) by some sequence of full turns of nearest 
neighboring threads if and only if the lattice sites corresponding to these threads belong to the 
same connected cluster of closed bonds. Otherwise these distant threads are not entangled. Thus, 
if a cluster of connected sites touches the opposite sides of the lattice, then we can always find two 
threads on opposite sites of the lattice entangled to each other. 

There are two formulations of the classical theory of bond percolation: canonical and grand 
canonical. In the canonical approach the averaging is performed over configurations with preserved 
number of closed bonds S on the lattice and, hence, the normalized concentration of closed bonds 
s = is also fixed. In the grand canonical approach only the probability p of a bond to be 
closed is fixed, while the concentration s fluctuates near the average value (s) = p. These two 
formulations are linked together in a simple way. If we know, for example, the crossing probability 
in the horizontal direction, TTfi{S), for the canonical distribution, we can obtain the same quantity 
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in the grand canonical formalism by averaging 7r/j(S') over the Binomial distribution 



N 



7rh{p,L) = Y,Pb{S;p,N)Trh{S,L) (2) 
5=0 

where Pb{S;p, N) is the probability to have exactly S closed bonds on the lattice of size L with 
N = 21? bonds for fixed value of the bond closure probability p: 

Pb{S;p,N) = (3) 

Here and below we assume that a function of the argument S (or of s) corresponds to the canonical 
distribution, while a function of the argument p corresponds to the grand canonical distribution. 
The capital letters S and M denote the number of closed bonds and number of dropped cells for 
the entire lattice, while the letters s = jj and m = ^ denote the same quantities normalized per 
number of bonds. 

In Refs. 3, 3) Q| the relation @ has been used for very efficient numerical simulation of the 
grand canonical crossing probability by using the canonical one, which is more fast and simple in 
practical applications. 

Let us note that eq.Q sets the general relation between two different ensembles (grand canon- 
ical and canonical) and we shall exploit this relation for our goals to establish connection between 
the standard grand canonical formulation of crossing probability and our canonical formulation of 
percolation produced by growing heap. Namely, in our model for each lattice size L and specified 
value r+ we randomly add each time moment the new cell to the heap and stop at the time moment 
when M cells are added. After that we check the horizontal spanning property of obtained config- 
uration of closed bonds by using the Hoshen-Kopelman algorithm J^]. Then we perform averaging 
over 10^ different configurations and obtain the topological crossing probability vrj(M; r+, L). The 
only difference from the usual crossing probability 7r/j(p, L) is that we generate each configuration 
of closed bonds by adding cells to the heap instead of closing each bond independently. The crucial 
difference appears for r+ 7^ 0; 1, i.e. when annihilation is allowed. Let us remind that for topo- 
logical percolation each closed bond is the projection of one (or several) cells from the heap to the 
base surface. 

In the FigEJwe have plotted the topological crossing probability TTt{m; , L) as a function of 
m = ^ for two cases r+ = 1; 0.5. In the same figure the crossing probability TTh{s;L) of the 
standard bond percolation is plotted for comparison as a function of s. The graph for the usual 
percolation 7r^(s; L) crosses the horizontal line 0.5 at the critical point Sc = \- The slop of vr^(s; L) 
is maximal at the critical point. The same construction can be used for the topological crossing 
probability. The graph iTt{m; r^, L) for r = 1 crosses the horizontal line Sc\ at the point with the 
maximal slop, and the graphs for different lattice sizes L cross each other at that point. As we 
see later, this is the critical point of topological percolation for r"*" = 1. The the similar behavior 
(however with different critical point) is valid for r = 0.5. 

We can express the topological probability TTt{M]r~^ = 1,L) via the crossing probability 
-Kh{M]L) in the way similar to eq. ((2]). Thus our next goal is to find the probability distribu- 
tion Pt{S;M,N) of having S closed bonds on the A^-bond lattice if we have added M cells to a 
heap. This probability distribution Pt{S] M, N) plays the same role for the topological percolation 
as the binomial distribution Pb(S;p, N) plays the role for the percolation in the "canonical" con- 
sideration. To begin with, let us consider only clockwise local turns (r"*" = 1) because in this case 
there are no cancellations of opposite turns, and the problem becomes essentially more simple. We 
denote this case "the irreversible topological percolation". 
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Consider some particular configuration on the lattice with S closed bonds and apply the mean- 



field consideration. Adding one extra cell to a heap, we see that with probability s 



N 



this cell hits 



closed already bond and does not change the number of closed bonds. However with probability 
1 — s the new cell hits the open bond and closes it increasing the number of closed bonds by one. 
Therefore, we can write down the master equation for Pt{S; M, N) on the lattice L: 



Pt{S+l;M + 1,N) 



^^Pt{S + 1; M, N) + (l-^] Pt{S; M, N) 



(4) 



with the following initial and boundary conditions: 

' p^(S = 0;M = 0,iV) = 1 
Pt(5 = 0;M > 1,A^) =0 
PtiS = l;M = l,N) = l 
Pt{S = I; M ^ 1, N) =0 

It is easy to solve this problem by the method of generating functions. Define 

oo oo 



(5) 



(6) 



Eqs. ©-© in terms of Q{t, z) read 

1 dQ{t, z) 
N dt 



M=Q 5=0 



1-tz 
tz{l - t) 



Qit,z) 



1 

tz{l - t) 



(7) 
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The solution of ((7| (properly normalized) is as follows 



N(i-z) IN N N 

Qit, z) = {l- t)^- 2F1 , 1 + AT _ 1 _ t 

z z z 



(8) 



where 2Fi{---) is the standard hyper geometric function — see, for example 20j. 



Now we can straightforwardly compute the mathematical expectation {S{M,N)) and the dis- 



persion {AS^{M,N)) = {S^{M,N)) - {S{M,N)y: 



(S(Af.iV)) = ^J^J^ 



dQ{t,z) 



{SiM,N)iSiM,N)-l)) 



1 r dz 



dt 

d^Q{t,z) 





-'1 




t=l. 







N ~ 1 
N 



M 



dt"^ 



t=l. 



M 



N I \ N 
where the contour C surrounds the point 2: = 0. Therefore 

{AS'^{M,N)) = {S{M,N){S{M,N) -1)) + {S{M,N)) - {SiM,N)f 



N{N -1) 



N -2 

N 



M 



N - 1 
N 



2M 



+ N 



N -1 
N 



M 



(9) 



(10) 



(S(M,N)) ^j^^ dispersion 



Let us evaluate the normalized expectation {s{m,L)) - 

(As^(m,L)) = ^^'^ ^'^^'^ , where N = 2L^, m = In the asymptotic regime, where 

^ 1, M ^ 1 and m = ^ = const, we can rewrite © for normalized quantities in the limit 
L — > 00 as follows 



{s{m,N)) 



and 



lim 



(5) 
N 



Nm 



1 — lim 

N—>oo 



N 



N' 



2Nm 



1 - e" 



Nm 



On the basis of ()12|) we expect for all r 



{As'^{m;r^,L)) < 



const const 



N 



2L2 



(11) 

(l + m)e-2™ 
(12) 

(13) 



Now we are in position to attack the "reversible topological percolation", i.e. the case r+ = 0.5 
when both clockwise and counterclockwise turns are allowed. First of all we estimate the maximal 
number of annihilated cells. In |3| the average number (H) of the most top cells (called "the 
roof) of the heap in the stationary regime M ^ N has been analytically computed. Only those 
cells are accessible for cancellations. We reproduce here the result of (13i] : 



{h{m,N)) 



{H{M,N)) 1 



N 



7 



(14) 
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For convenience we have reproduced in Appendix the derivation of this expression. 

Following the methods developed in the works [ll,!!^,!!!] one can also compute the dependence 
{h{M,N)) in the intermediate regime M < N. The function {h{m,N)) in all regimes is well 
approximated by the curve 



{h{m,N)) = -{l 



-7m\ 



(15) 



The analytic derivation of this expression is presented in Appendix. To verify (|15j) we compare in 
Figinithe data of numerical simulations for the average size of the roof {h{m; r+, L)) for r+ = 0.5, 
L = 128 with eq.ldSJ). One sees that the function {h{M, N)) = \{1 - e"^'") perfectly fits the data. 
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Figure 6: Dependence the size of the roof of the heap h{m; r+, L) for r+ = 0.5 and L = 128 in comparison 
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Let us make one step back and demonstrate that for r"*" = 1 we can arrive at eq. (|11() in a simple 
way. Suppose that some configuration of cells covers the base surface. If we put one more cell, it 
hits the empty bond and increases the number of closed bonds by one with the probability 1 — s, 
while it hits the closed bond and does not change the number of closed bonds with the probability 
s. These conditions are summarized in the table below: 

{+1 with the probability 1 — s 
(16) 
with the probability s 

So, we have 

' d{s{m,r = 1)) ^ . ^ 

; = 1 — Sim) 

dm (17) 

s{m = 0) = 
The solution to this equation is ((TT|l . 
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We modify now eq. ()17|) to catch the "reversible topological percolation" (r"*" = 0.5), taking into 
account the annihilation. We have stressed that the annihilation of cells occurs in the roof only. 
Each bond of the roof is either black or white with the probability ^. Hence the total probability 
of annihilation is about ^ {h{'m; , L)) . Moreover, the bond in the base surface becomes "open" 
if and only if the cell of the roof is single in the corresponding column, i.e. below the given roof's 
cell there are no other cells in the same column. The probability to have empty set of 1 — (s(m)) 
columns is given by the solution of ea. (|17|) : 1 — {s{m)) = e"*". Hence, the probability Pcond of 
annihilation of the roof's cell under the condition that below this roof's cell there are no other cells 
in the same column (assuming the uniform distribution of the roof's cells and empty columns), is 
given by the product: 



Pcond = 7j {Hm; r+ = 0.5)) x 



Therefore 



14 



-(1-e 



-7m\ 



(18) 



As(m, = 0.5) = < 



-7m\ 



—1 with the probability %^(1 — e~ 

with the probability s — ^xj-(l — e~'^™') 
+1 with the probability 1 — s 



(19) 



Hence, we can write 



f d(g(m),r+=0.5) / ^ , 

= 1 — sm 1 — e 

dm ^ ^ 14 ^ ' 

s{m = 0) = 



(20) 



The solution to eq. (^0]) reads 



(s(m;r+ = 0.5)) = 1 



-8m 



97 e""" me- 



98 98 
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(21) 



In FiglTJwe plot the numerical data for {s{m;r^ , L)) for = 1,0.5 and L = 128 (crosses). We 
see, that data for r+ = 1 is in the excellent agreement with eq. (|11() . For comparison, in the same 
figure we plot also the analytic expression 1)21(1 for the "reversible topological percolation". As one 
sees, this curve fits rather well the corresponding numerical simulations. 

Let us introduce the probability distribution for normalized quantities Pt{m,s,L) = 
2L^ Pj(2L^s, 2L^m, 2L^). We can write (at least for the critical region) the scaling expression 
for the crossing probability 

7r,(.,L)^/((.-s,)Li/'^) (22) 
where = | is the critical exponent of the correlation length for the percolation. Therefore 

c^L'/'^f"{{s-s,)L'/n (23) 
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Figure 7: Comparison of (s{m] r+, L)) for r = 1, 0.5 and L — 128 with the approximation functions. 



For 5 » 1, iV » 1 we have 



N 1 

(m„L) = y2Pt{S;Mc,N)7rh{S/{2L^),L) ^ / Pt{s;mc, L)Trhis, L)ds 
5=0 -^0 



dTTh{s,L) 



ds 



As + 



As^ + o(As2) ds 



TTh{Sc,L) + 



d'^TTh{s,L) 



ds"^ 



{As\m, L)) ^ TThis. L) + L^'^-^f"{s - s,)) + o{L^I^~^) 

'c 

(24) 

where As = s — Sc, Sc = 1 — e "^^ and the function /(s — Sc) is regular in the vicinity of the point 
Sc- We use the condition H13|) for the dispersion (As'^{m,L)\. Hence for large lattices one has 



lim TTt{m; r = 1, L) = vr/i ( 1 — e L 

L^oo 



(25) 



The bond percolation on infinite square lattice occurs at the concentration of bonds Sc = Pc = ^■ 
Hence, using ((TT|). we get the critical concentration of local turns (or falling cells) rric = ^ at the 
percolation threshold in the limit L ^ 1 for irreversible ((r"*" = 1) and reversible (r"*" = 0.5) cases: 

mc(r+ = 1) = ln2 ~ 0.693 for r+ = 1 

(26) 

mc(r+ = 0.5) « 0.734 for r+ = 0.5 



These solutions are obtained from eqs. (|11() and H21|) correspondingly: 
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Figure 8: The topological crossing probability 7rj(s; r+, L) as a function of the average number of the closed 
bonds s — I — exp(— m) for r = 0.5, 1 and crossing probability 7r/i(s; L) of bond percolation, canonical case. 



The values (|'2f)j) are in good agreement with the data of numerical simulations shown in FigEl 

In the FiglSl we show by symbols the topological crossing probability vr^ as a function of a 
variable s = 1 — exp(— m) for = 1,0.5 and L = 32,64,128. In the same figure we plot the 
crossing probability of bond percolation TTfi{s; L) for the canonical case by lines. One sees that the 
symbols for r+ = 1 lie on the appropriate lines. Our investigation allows us to conclude that the 
topological percolation belongs to the universality class of the two dimensional percolation. We 
confirm this statement numerically. Namely, we extract the correlation length exponent = 4/3 
from the numerical data shown in the Fig|5J The outline of our construction is as follows. The 
eq. (|22|) allows us to conclude that the derivative of the crossing probability at the critical point 
scales with the lattice size as 

<(.„L)^/'(0)lV'^ (27) 

Now we define numerically the derivative of the crossing probabilities. We approximate the data 
for 7rt(m; r"*", L) in the vicinity of the critical point 0.5 — 0.05 < Trt{m; r"*", L) < 0.5 + 0.05 by the 
linear function 0.5 + a(r, L){m — mc(r+, L)). We perform this procedure for grand canonical and 
canonical distributions of usual percolation (crossing probabilities TThip, L) and '7r/i(s, L)) as well as 
for the topological percolation = 1,0.5 (crossing probability Tit). Then we plot the derivative 
of the crossing probability at the critical point TT'{sc,r,L) = a{r,L) as a function of the lattice 
size L for grand canonical and canonical distributions for the standard percolation, as well as for 
irreversible (r+ = 1) and reversible (r+ = 0.5) topological percolation. The corresponding results 
are shown in FiglHl 

The data depicted in Figl^lare approximated by the function aL^/^ . The result of approximation 
is plotted in the FigElby lines. One sees that the numerical values of the inverse correlation length 
length 1/v coincide with the analytical value | within accuracy of the approximation. 
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Figure 9: Extraction of the inverse critical exponent of the correlation length form the derivative of the 
crossing probability at the critical point. 



IV. CONCLUSION 



We have investigated the topological phase transition in the bunch of randomly entangled 
directed threads. Specifically we pay the most attention to the determination of the minimal 
number of local turns necessary to produce the fully topologically connected cluster of threads. 
Namely, below Mc the opposite side of the square lattice of threads are disjoined, while above 
the opposite sides belong to the single cluster of connected threads. Two models of topological 
percolation are considered: "irreversible" and "reversible". In the irreversible case (r"*" = 1) all 
local turns are only clockwise, while in the reversible case both clockwise and counterclockwise 
local turns are available with equal probabilities {r~^ = 0.5). 

We map the problem of topological percolation onto the standard two-dimensional percolation 
and relate the above defined value Mc (normalized per the number of lattice bonds) to the critical 
value pc of the percolation threshold on the square lattice. The consideration of the reversible 
topological percolation demands special care. Wc give the corresponding estimate for the critical 
value Mc considering the reversible topological percolation as the growth of the random heap of 
pieces. In particular we estimate the probability to "open" the bond of the base surface from the 
probability of cancellation of a piece in a growing heap. 

In addition, we find numerically the critical exponent for the topological percolation and show 
that is is equal the correlation length exponent u = ^ of the standard two-dimensional bond 
percolation. 



Appendix A 



The process of growth of a heap (i.e. the random walk on the surface locally free group CT^d) 
consists in adding step-by-step new "black" or "white" blocks to the roof. The dynamics of a 
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heap is controlled by the dynamics of a roof. For a particular configuration of a heap we define 
the number H of most top segments in a heap (i.e. the "size of a roof) as well as the number 
"non-roof's" segments, Li, having exactly i neighboring roof's segments — see (l^ . (Remind that 
there are N = 2L^ lattice bonds on the lattice). For the values H, Li the following conditions hold: 



\ L0 + L1+L2 + H = 2L'^ + 2L 
[dH - 8(L + 1) < Li + 2L2 < 6H 

If L ^ 1 one can neglect the boundary conditions and rewrite (|A1|) is simpler form 

( L0 + L1 + L2 + H = 2L^ 
\ Li + 2L2 = GH 

For a given configuration the local dynamics of a size of a roof reads 

AH = 1 with the probability ^ 
AH = with the probability ^^^^ 
AH = — 1 with the probability 

Equations ()A3|) allow to write the following equation for the expectation {H{M)): 

d{H{M)) L0-L2 



(Al) 



(A2) 



(A3) 



dM 



2L2 



This equation can be rewritten for the normalized quantities {h{m)) = i^ji ^ind 
closed form with the help of ()A2|1 . We get: 

f d{h{m)) 



m 



M 
2l7 



1 - 7h 



The solution to HA4() reads 



dm 

{h{m = 0)) = 
1 



in the 



(A4) 



{h{m)) = -(1 



-7m\ 



(A5) 



The derivation of ()A5I) is the basis for the approximation H15|) . In a stationary case (i.e. for 
m — > 00) we have {h(m)) = i. This expression coincides with dJ. 
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